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The evolution of piecewise constant distributions of a conserved quantity related to the frozen-in 
canonical vorticity in effectively two-dimensional incompressible ideal EMHD flows is analytically 
investigated by the Hamiltonian method. The study includes the case of axisymmetric flows with 
zero azimuthal velocity component and also the case of flows with the helical symmetry of vortex 
lines. For sufficiently large size of such a patch of the conserved quantity, a local approximation 
in the dynamics of the patch boundary is suggested, based on the possibility to represent the total 
energy as the sum of area and boundary terms. Only the boundary energy produces deformation of 
the shape with time. Stationary moving configurations are described. 



PACS numbers: 52.30.Cv 

I. INTRODUCTION 

In this paper we consider a special class of vortical flows 
in plasma that correspond to the model of ideal electron 
magnetohydrodynamics (EMHD) (see, e.g., |], @, §, £| 
and references therein about EMHD and its applica- 
tions). In general, we deal with ideal EMHD-flows that 
are effectively two-dimensional (2D) and completely de- 
termined by a single function of two spatial coordinates 
and the time, namely with the usual planar flows, with 
the axisymmetric flows having the only azimuthal compo- 
nent of the magnetic field, and with the flows possessing 
the helical symmetry of frozen-in lines of the generalized 
vorticity. Our purpose here is to introduce a simplified 
analytical description adopted for the case when the func- 
tion is piecewise constant (such distributions are usually 
called "patches"), with a typical spatial scale between 
well separated the electron inertial length and the ion in- 
ertial length. In a narrow region near sharp boundary of 
a "patch" the density of the electric current is high, thus 
the boundary can be considered as a dissipationless cur- 
rent sheet. From analytical viewpoint, patches are the 
simplest examples of EMHD-flows with current sheets, 
in the sense that no other degrees of freedom are present 
in the system, except the shape of the patch. It is a well 
known fact that spontaneous formation of current sheets 
from initially smooth magnetic field configuration is a 
typical phenomenon in plasma dynamics, and the role of 
current sheets is very important. That is why it is inter- 
esting to investigate the problem of their motion, even in 
such highly idealized formulation as we do in this work. 
As an explicit instance of fundamental physical phenom- 
ena tightly linked with our analytical study, we can in- 
dicate the fast magnetic field penetration into plasmas 
due to the Hall effect, since in some conditions the ideal 



EMHD is applicable to describe this phenomenon j|, [|. 
More precisely, we focus on the effect of the electron in- 
ertia on the dynamics of a narrow moving front between 
magnetized and unmagnetized plasma for those special 
highly symmetric configurations of the generalized vor- 
ticity field. It should be noted that on the mentioned 
scales the magnetic field is similar to the generalized vor- 
ticity field almost everywhere except the close vicinity of 
the sharp patch boundary, where the magnetic field is 
smoothed on the electron inertial length. 

As known, the EMHD model approximately describes 
the motion of the low-inertial electron component of 
plasma on sufficiently short scales (below the ion iner- 
tial length), while the much heavier ion component may 
be considered as motionless |jT| . Conditions for electrical 
quasi-neutrality of the plasma are assumed, which imply 
that the macroscopic velocity v(r, t) of the electron flow 
should be sufficiently small and also all potential waves 
and oscillations of the plasma density should not be ex- 
cited. Only vortical degrees of freedom of the system 
are relevant in such circumstances. With these condi- 
tions, everywhere in the space the concentration n(r,t) 
of electrons is approximately equal to the prescribed con- 
centration N(r) of ions, thus the density of flow of the 
electrons, 

j(r,i) = ra(r,t)v(r,i), 

is almost divergence free, divj » 0. The role of ions 
is thus reduced simply to compensation of the electric 
charge of electrons by providing the static neutralizing 
background. Such flows of the electron fluid create the 
divergence-free field of electric current density — ej, where 
— e is the electron charge. Hence, the quasi-stationary 
magnetic field 

B(r,t) = -— curr^r.t) (1) 
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should be taken into account when considering forces act- 
ing upon the electron fluid. 
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Important for our consideration is that in some cases 
dissipative effects due to finite resistivity and/or viscos- 
ity may be neglected without dramatic loss of accuracy. 
So we study below the conservative EMHD model. A re- 
markable feature of the conservative hydrodynamic-type 
systems is that all of them possess an infinite number of 
integrals of motion related to a basic physical property 
of fluids, the relabeling symmetry (see, e.g., @, ||, [nj 
and references therein for a recent discussion) . This sym- 
metry manifests itself as the freezing-in property of a 
divergence-free vector field (the canonical vorticity field) 



n(r, t) = curl 



6C 



(2) 



specified by the Lagrangian functional £{n, j} of a model 
and evolving with time accordingly to the equation 



ft t = curl[v x n]. 



(3) 



For example, in the usual Eulerian hydrodynamics J~i(r, t) 
is simply proportional to the velocity curl. As to the ideal 
EMHD, it is well known fact and we will see it below once 
more that frozen-in is the field 

f2(r, t) = curlmvfr, t) B(r, t) 

c 

Aire 2 i. 
= curlmvH — curl j, (41 

cr 

where m is the electron mass. 

Below we deal with a spatially homogeneous ion back- 
ground N = const, thus the flows of the electron fluid are 
supposed to be incompressible, (V • v) = 0. There ex- 
ist three classes (hereafter referred as P, A, H) of incom- 
pressible ideal EMHD flows having a special symmetry 
which allows one to describe the vorticity distribution in 
terms of a single scalar function uj depending on two spa- 
tial coordinates and the time. Besides the usual planar 
flows with 



ftp = e z u>(x,y,t), 



(5) 



these effectively two-dimensional are the axisymmetric 
flows with zero azimuthal velocity component, 



where 



& A = [e z x r]uj(z,q,t), 



q=(x 2 +y 2 )/2 



(6) 



and also the flows with the helical symmetry of the 
frozen-in vortex lines, 

fl z H — uj(xcosKz + ysmKz, —xsinKz + ycosKz,t), 

(7) 



fl x H = -KyQ z H , n y H = Kxfl z H , 



(8) 



that are space-periodic along z-direction with the period 
L z = 2%/K. What is important, in each of these three 



cases the evolution of the corresponding ui(u,v,t) is noth- 
ing else but only the transport of its level contours by 2D 
incompressible flows. The general structure of the equa- 
tions of motion for the function u)(u, v,t) is 



0, 



(9) 



with the stream-function ^(u, v, t) being specified by the 
Hamiltonian functional H. a {uj(u, v)} of the model, de- 
pending also on the type of symmetry (a = P, A, H), 



* = 



Soj 



(10) 



We will see later that the Hamiltonians Ti a {uj(u, v)}, cor- 
responding to the ideal EMHD, take the quadratic form 

7-C a {u)(u, v)} = — I uj(u,v)G a ui(u,v)du dv, (11) 



where the nonlocal linear operators G a possess smooth- 
ing properties like the usual two-dimensional A -1 - 
operator. Therefore the flows with discontinuous piece- 
wise constant distributions of the function uj(u, v) are 
possible. From Eqs.(PHl0|) it follows that the shape 
t), u(£, i)} of boundary of such a patch with u> = 
luq = const (here £ is an arbitrary longitudinal parame- 
ter) , evolves in accordance with the equation 



U t V£ - U£V t 



l_d_ 

uj <9£ 



(12) 

where {«(£), w(£)} is the energy of the patch expressed 
through its shape. The purpose of this work is to study 
the motion of the patches in (u, w)-plane for all three 
kinds of the geometric symmetry. Wee will extensively 
use the fact that in EMHD the kernels G a {u\, V\\ U2, 1*2) 
of the operators are exponentially small if the distances 
between (ui,v±) and (1*2, V2) are much longer than an in- 
ternal length parameter of the problem, the inertial elec- 
tron skin depth d = (mc 2 /47re 2 n) 1/ ' 2 . For comparatively 
large patches, this property makes possible to single out 
the local boundary term in the Hamiltonian, which is re- 
sponsible for the evolution of the patch shape, while the 
main area term in the approximate local Hamiltonian 
only results in a uniform motion of the patch without 
changing the shape. 



II. CANONICAL FORMALISM FOR EMHD 

For convenience and self-consistency, now we repro- 
duce briefly derivation of the EMHD equations, following 
the canonical formalism adopted for fluids as described 
in Refs. [|, [l(| As the start point, let us consider the mi- 
croscopic Lagrangian of a system of electrically charged 
point particles that can be written up to the second order 
on v/c, as it is given in the famous book by Landau and 
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Lifshitz [hi 
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-y 

1 „2 



e a eb 



nucro / . 2 2 _ |ra _ r 



^ 8c 2 



e a e 6 



4c 2 ^-f |r a - r b 



{v Q • v b + (v a • n Qb )(v b • n Qb )](^3) 



where r a (i) are the positions of the point charges e a , 
v a {t) = r a (t) are their velocities, n a b(t) are the unit vec- 
tors in the direction between e and e b . The first double 
sum in Eq.(|l3"|) corresponds to the electrostatic interac- 
tion, while the second double sum describes the magnetic 
interaction via quasi-stationary magnetic field. It is very 
important that for a system with macroscopically large 
number of particles the magnetic energy can be of the 
same order (or even larger) as the kinetic energy given 
by the first ordinary sum in Eq.([l3|), while the terms 
of the fourth order on the velocities are often negligi- 
ble. Generally speaking, a large part of plasma physics is 
governed by this Lagrangian, at least in the cases when 
the velocities of particles are non-relativistic and the free 
electro-magnetic field is not excited significantly. Obvi- 
ously, different physical problems need different proce- 
dures of macroscopic consideration of this system. The 
most accurate (and the most complicated) would be a 
kinetic description. However, for our purposes it is suffi- 
cient to apply more simple and naive procedure of the hy- 
drodynamical averaging, that gives less accurate descrip- 
tion of the system in terms of the concentration n(r, t) of 
electrons and the density j(r, i) of their flow, that satisfy 
the continuity equation 



on .. . 
- + d lvJ = 0. 



(14) 



Neglecting all dissipative processes that take place due to 
collisions of the particles (though on this step we strongly 
reduce applicability of the following conservative EMHD 
model) , and considering the ions as macroscopically mo- 
tionless, we derive from Eq.(|l3|) the following Lagrangian 
functional of the electron fluid on the given static ion 
background with the macroscopic concentration iV(r): 



£ e i>J} = dr 



m- - s(n) 



dridr 2 



e 

~4? 



ri - r 2 | 
dri<ir 2 



■{n( ri ) - N( ri )}{(n(r 2 ) - N(r 2 )} 



j( r i) -j( r 2) 



|ri - r 2 

(j(ri) • {ri - r 2 })(j(r 2 ) • {n - r 2 }) 



ri - r 2 



,(15) 



where the internal energy e(n) of the electron fluid takes 
into account energy of the thermal disordered motion 
of electrons and also the microscopic-scale-concentrated 
part of the electro-magnetic energy. It should be kept in 
mind that s{n) depends also on the specific entropy but 
we suppose the flows isentropic. 



Since we are interested below in relatively slow vortical 
motion of the electron fluid, when all possible potential 
waves and oscillations are not excited significantly, so the 
quasi-neutrality condition n ss N is well satisfied, it is 
possible to neglect the second line in the expression JlE|), 
as well as variation of e(n). Thus, for the EMHD, which 
describes this special dynamical limit of the system ( |l5| ) , 
we have the Lagrangian 



£n{j} = 



Aire 2 (curl" 1 ]) 5 



2n 



dr, 



(16) 



where the concentration n(r) is a prescribed function, so 
the density j of flow is divergence-free, divj = 0. 

In general, the equation of motion for the divergence- 
free flux field j has the form (compare with |l(| ) 



di cml {-IF ) = curl 



- x curl — -— 
n \ Sj J 



(17) 



In the particular case for EMHD this yields 



d 



mj Aire* 



— curl — + 
at \ n 



curl 



curl 



curl — 

n 



47re" 



curl j 



One can consider the canonical vorticity field 

Si 



$7 = curl 

and the Hamiltonian functional 
' 5C n 



Using the relation 



Si 



j dr - £„ 



j=j{":«} 



j = curl 



5H 



.(18) 



(19) 



(20) 



(21) 



it is possible to rewrite the equation of motion ( |l7| ) in 
the form 



fit — curl 



nirl 1 m) x S 



(22) 



that emphasizes the freezing-in property of the canonical 
vorticity. 

Hereafter we consider the case n — const and use di- 
mensionless variables, with all length scales normalized 
to the ion inertial length. The expression for the EMHD 
Hamiltonian can be written then in the following form 

H{n} = ^ J J G(|r 1 -r 2 |)(0(r 1 ).fi(r 2 ))rfr 1 dr 2 , (23) 



with the Green's function 



G(r) 



exp(— rj A) 
47rA 2 r 



(24) 
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FIG. 1: 
flows. 



Sketch of generalized vorticity patch in axisymmetric 



where A = \J m/M is a small parameter, the ratio of 
the electron inertial length to the ion inertial length. It 
should be noted that the EMHD model in this simplest 
form, without taking into account the ion motion, is ap- 
plicable only on scales below the ion inertial length. A 
more accurate model differs by another Green's function 
in the double integral (H|) : instead of the expression ( |24| ) 
with effectively finite radius of interaction, one has to use 
the modified Green's function, with infinite radius of in- 
teraction (see, e.g., [[L2] for more detail) 



G(r) 



_1_ / cxp(-r/A) | 1 
47r V \ 2 r r 



(25) 



However, in this work we deal with the function (24) 



since our goal is to construct a local approximate model. 
Thus, a typical size L of vortex structures must be in 
the limits A <C L <C 1, in order the contribution to the 
Hamiltonian from the second term in Eq.(p5|) to be small 
in comparison with the contribution from the first one. 



III. LONG-SCALE LOCAL APPROXIMATION 
IN CONTOUR DYNAMICS 

A. Axisymmetric flows 

Let us consider effectively 2D flows, for instance, the 
axisymmetric flows (^). The dynamics of the function 
u)(z, q, t) is determined by Eqs.(||- [!(]), where TLa{^} — 
(l/2n)'Hiuj\e z x r]}. For long-scale flows, the Green's 
function (|24| ) in the double integral (|2^) is almost the 
same as the (5-function. Therefore, large vortex patches 
like that shown in the Fig.[j], with constant u> and sharp 
boundaries, can be approximately described by a local 
Hamiltonian broken onto two parts — the bulk energy 



njD 1 

riA = - 



n rdr dz = i 



D 



, Q i qdqdz, 

' D 



and a "surface" energy originated by the effect of non- 
locality near the boundary. In leading order the bound- 
ary term can be calculated as if we have locally a ID 
configuration with the jump Q = ujQ^/2q in the vorticity 
field. In such a ID case, the additional energy per unit 
area of the boundary is simply 



n 2 

2- 2 



-CA 



dC 



r> 2 A 



that gives the surface energy of the patch in the axisym- 
metric flows, 



\UJn 



H% D = i r 2 ■ r^/(dz) 2 + (dry 

4 J 3D 



1 J 3D 



z) 2 + [dqf 



(26) 



Let us for simplicity take u>o = 1. Using the explicit ex- 
pression for the Hamiltonian functional H* A {z(t;), q(£)} = 
H° + H d A D , 

WaMO, 1(0} = j> qzq^ ~^ j> Q\J^QZ 2 + q 2 d^, 

(27) 

we obtain from Eq. (H2) the following equation of motion, 
which does not depend on choice of the longitudinal pa- 
rameter £: 



z t q£ - Ztft = 95 + A 




(28) 



This local nonlinear equation is one of the main results 
of present work. It approximates the nonlocal contour 
dynamics in axisymmetric EMHD-flows if a typical lon- 
gitudinal scale L (in (r, z)-plane) of the contour satisfies 
the condition A<L< 1. 

Now we study stationary moving (along z-axis) con- 
figurations. Let the shape of a patch loq — 1 in axisym- 
metric flow be given (locally) by the function z(q, t) (this 
implies the fixed choice of the longitudinal parameter in 
Eq. (p8[) , <? = £). Then the corresponding equation of 
motion is 



z t (q) = I + A 




(29) 



Let us consider stationary moving solutions, z t = 1 + 
2AC, where C = const. The shape of the boundary in 
such case is determined by the ordinary differential equa- 
tion 




= 2C, 



that can be easily integrated: 



z = I dq 



[Cq 2 + dq + b) 



2q(Cq 2 + dq + bf 



(30) 



(31) 
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FIG. 2: The parametric plane for axisymmetric stationary FIG. 3: Examples of axisymmetric solutions: crimped current 
solutions. channels and a magnetic ring. 



with some constants a, 6, C. To simplify further analysis 
it is convenient to make the rescaling 



aq, 



az. 



1 

2C72' 



In these new coordinates the stationary solutions depend 
only on two parameters, 



aq + b 



q{q 2 + aq + bf 



(32) 



However, in general case the constants a and b may 
not take arbitrary values, since the expression under the 
square root in the denominator must be positive in some 
range of (positive) q. Another restriction is that self- 
intersections of the curves are forbidden. The parametric 
(a, 6)-plane is shown in Fig.|[ First, there exist z-periodic 
configurations that occupy 2D region in the parametric 
plane and correspond to axisymmetric current channels 
with a sharp crimped surface. On different pieces of the 
boundary of that region we have different types of sta- 
tionary solutions. Partially the boundary is determined 
by the parametrically given curve (the thick solid line in 
Fig.§ 

o(r) = -2r(r-3/4), b(r) = r 3 (r - 1/2). (33) 

If the parameter r is within the limits 3/8 < t < +oo, 
then the solutions near the curve look like linear sinu- 
soidal waves with vanishing amplitude. If < r < 3/8, 
then z-period of the current channels tends to infinity, 
and we have highly nonlinear solitary waves. For the 
third piece of the boundary, where self-touch of the cur- 
rent channels with multi-valued dependences q{z) takes 
place, we do not have a simple analytical expression, how- 
ever it can be found numerically. Second, besides the 
infinitely long current channels, there exist finite con- 
figurations (cross-sections of magnetic rings). For finite 
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FIG. 4: Cross-sections of magnetic rings. 



configurations the curves in (z, q)-plane must be closed, 
thus in this case the integral ( p2[ ) between the two cor- 
responding zeros of the denominator must be equal to 
zero, that gives a relation between a and b (the numer- 
ically calculated separate curve in Fig.|). Though the 
piece of the parametric curve (|3^) with — oo < r < 
corresponds to nothing in the parametric plane, but the 
line of finite configurations approaches it from below, as 
a -> 0. 

Examples of appropriate solutions are presented in 
Fig.|| and Fig.[I| Strictly speaking, only the magnetic 
rings can satisfy the condition L« 1. For the infinitely 
long current channels, applicability of the local model is 
not so well justified as for the rings. 

In the limit, when a — > 0, the "upper" part of the con- 
tour shape is the parabola when being viewed in (q, z)- 
plane (the half-circle in (r, z)-plane), 



(z 2 + r 2 ) 
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For the rings, the tendency towards this shape is clearly 
seen in the FigJJ, as \a\ decreases. 

B. Helical flows 

Analogously, a patch u = 1 in the helical flows can be 
investigated. For this case the bulk energy is 

2 
1 



njD 



[1 + K 2 (u 2 



J )]du dv 



D 



J R2 I < 2 

udv H (p (uv 

2 JdD 2 / aD 



u 6 /3)dv, (34) 



[l+X 2 w 2 ]^/(dw) 2 + K 2 (wdw) 2 (35) 



while the boundary energy can be expressed as follows: 

njdD ^ 

Here the square root is proportional to the area element 
of the helical tube surface. 

Note that the planar flows (^]) are included into this 
consideration, since they correspond to the case K = 0. 
The boundary term in that case is simply proportional 
to the length of the contour that results in the integrable 
dynamics fll^ , |l4f| — the evolution of the boundary cur- 
vature k(1, t), where I is the arc-length along the contour, 
is determined by the modified KdV equation, 



Km 



-K 2 Ki 



(36) 



In the case K ^ 0, the term just produces the 
uniform rotation of the contour with the angular ve- 
locity (fio — —K 2 around the origin. The shape evo- 
lution generated by Ti.jp is more convenient for ana- 
lytical study in the polar coordinates (note the relation 
d{u,v)/d(q,ip) = 1) 

u = \f2q cos <p, v — y^2q sin tp, (37) 

since the expression for TLjf{ip(q)} does not contain <p(q) 
itself, but only d q <p(q): 



K d H D = ~ I dq[l 



A 



2K 2 q]^\/2q- 
The equation of motion is 

dq \ 



2<7^ 2 



K 2 (38) 



-K z -- 




2q^ q [l + 2K 2 q] 



2qtf 



'l/2q , _, r(/ 

For stationary rotation, when ip t — —K 2 — 

(Cq 2 + Bq + A)y/l + 2K 2 q 



(39) 



+ K 2 

CA,we have 



dq 



Q ^/2q{l + 2K 2 q) 2 - A{Cq 2 + Bq + A) 2 

(40) 

Here also a right choice of the constants A, B, C must 
satisfy the condition for the curve to be closed and do 
not have self- intersections (see Fig.[|). 




FIG. 5: Examples of stationary rotating patches in helical 
flows. 



IV. DISCUSSION 

Thus, a local approximation for a class of effectively 
2D flows in the ideal EMHD, based on the Hamiltonian 
formulation of the problem, has been suggested for the 
nonlinear contour dynamics. Stationary moving config- 
urations with current sheets have been analytically (in 
integral form) found within this approximation. They 
describe propagation of magnetic rings, traveling waves 
on a crimped surface of current channels, including soli- 
tary waves, and rotating helical magnetic structures. Of 
course, our results are not able to explain all the impor- 
tant phenomena in EMHD. A payment for the relative 
simplicity of our local models has been the general im- 
possibility to consider in such a manner short-scale (~ A) 
contour perturbations, though their dynamics may have 
influence on the problem of stability of the stationary 
moving configurations. Also, the question about evolu- 
tion of 3D perturbations of the magnetic rings, current 
channels, and/or helical structures cannot be answered 
in such very simplified description of the flows by mov- 
ing curves in a plane. But we believe our approach useful 
since it develops an insight in qualitative theoretical un- 
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derstanding of such complicated nonlocal theory as the 
EMHD. For instance, we do not see another simple and 
compact way how to derive equations of motion for cur- 
rent sheets, analogous to Eqs.(p8|) and (|3^). 
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